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ABSTRACT 



We report on the discovery of a galaxy cluster ai z — 1.490 originally selected as an extended X-ray source in the 
'KMM-Newton Distant Cluster Project. Further observations carried out with the VLT-FORS2 spectrograph allowed the 
spectroscopic confirmation of seven secure cluster members, providing a median system redshift of z = 1.490±0.009. The 
color magnitude diagram of XMMU J0338. 8-1-0021 reveals the presence of a well populated red sequence with z— H « 3, 
albeit with an apparent significant scatter in color. Since we do not detect indications for strong star formation activity in 
any of the objects, the color spread could indicate different stellar ages of the member galaxies. In addition, we found the 
brightest cluster galaxy in a very active dynamical state, with an interacting, merging companion located at a physical 
projected distance of d ~ 20kpc. From the X-ray luminosity we estimate a cluster mass of M2oo~1.2 X 10^"* Mq. The 
data seem to suggest a scenario in which XMMU J0338.8-f 0021 is a young system, possibly caught in a moment of 
active ongoing mass assembly. 

Key words, galaxies: clusters: general - X-rays: general 



1. Introduction 

The search for distant (z > 0.8) galaxy clusters has ex- 
perienced a remarkable boost recently. By means of var- 
ious selection methods, several clu sters have been found 
in t he redshift de s ert, a t z > 1.5. iPapovich et al.l (|2010[ ) 
and iTanaka et al.l (j2010[ ) independently reported on a new 
IRAC selected cluster at z = 1.62, characterized by a well 
populated and tight red sequence. 

Further multi-band analysis revealed that the star for- 
mation activity in this cluster increases with the environ- 
ment density, with star forming gal axies being main ly lo- 
calized in the center of the system (jTran et al.ll2010[ ). An 
analogous trend has been observed in three other X-ray se- 
lected galaxy clusters, two newly confirmed a t z = 1.56 
(|Fassbender et al.ll201ll) and z = 1.58 ( Santos et al.ll201lf ) 
and one at z = 1.46 (|Hilton et al.l [20091 |2010j). For these 



* Based on observations under programme ID 084.A-0844 col- 
lected at the European Organisation for Astronomical Research 
in the Southern Hemisphere, Chile, and observations col- 
lected at the Centro Astronomico Hispano Aleman (CAHA) 
at Calar Alto, operated jointly by the Max-Planck Institut 
fiir Astronomic and the Instituto de Astrofi'sica de Andalucia 
(CSIC). 



systems, however, the galaxy red sequence is not fully es- 
tablished yet, revealing that the members are still building 
their stellar mass via star formation activity. The situation 
is different in clusters found so far at z < 1.4, where ac- 
tive star forming galaxies preferential l y reside in the clus- 
ter outskirt s (e.g . iBaloeh erall 120041: iLidman et all 120081: 
iPatel et al.ll2009t) . However, it is expected that more dis- 
tant galaxy clusters exhibit a larger n umber of star forming 
galaxies in their ce ntral regions (e.g. iHopkins et all 120081: 
iRettura et al.|[2010| ). In addition, early- type galaxies (ETG) 
in clusters are known to have completed the formation of 
their stellar populations much faster with respect to ETG 
in the field and the red sequence is expected to fade by 
z w 2 (jCobat et al.ll2008[) . The latter predic tion IS consis- 
tent with the recent results in iGobat et al.l (j201lh . where 
the discovery of the most distant cluster known to date, 
at z = 2.07, is reported. In this system the population of 
passive galaxies shows a very large intrinsic scatter in the 
Y — Ks colour, with an apparent absence of an established 
red sequence. This suggests that observations are now ap- 
proaching the early stages of galaxy formation in the dens- 
est environments. 

All these new findings highlight the importance to de- 
tect new distant galaxy cluster study targets for a syste- 
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Fig. 1. Left: H-band image (4.5' side length) of the environment of the galaxy cluster XMMU J0338. 8+0021 aX z ~ 
1.490. The XMM-A^ewton detected X-ray emission is shown by the log-spaced blue contours, with the four lowest ones 
corresponding to significance levels of 2, 2.7, 3.6, 4.6 a above the mean background. The associated ovcrdcnsity of color- 
selected red galaxies is displayed by the red contours, showing linearly spaced significance levels of 2-20 a. The big, black 
circles indicate the 0.5'/l' radii around the X-ray center; the small black ones mark the spectroscopically confirmed 
cluster members, with the corresponding IDs. The image was smoothed with a 3 pixel (1.3") Gaussian kernel. Right: 
z-|-H-band color composite image of the core region of the cluster. X-ray contours are displayed in cyan, red galaxies 
with 2.6 <z— H< 3.5 are marked by red dashed circles, spectroscopic members by white circles and interlopers (i.e. 
spectroscopically confirmed non-members) by green crosses. The 30" radius from the X-ray center is indicated by the 
yellow circle. The white ruler represents the beam size (FWHM) oi^MM- Newton at the observed off-axis angle of ^5.6'. 



matic investigation of those processes which drive galaxy 
evolution in different environments. 

In this Letter we report on the discovery of an X- 
ray selected, spectroscopically confirmed galaxy cluster at 
z=1.490, found within the XMM-Newton Distant Cluster 
Project (XDCP). We first present the X-ray data (glT|) 
followed by a discussion of the near-infrared (NIR) imag- 
ing observations ( ^12.21) and the optical spectroscopy results 
(©. In Section|3]wc will discuss the dynamical state of the 
system and its galaxy population. Finally, a summary and 
conclusions are reported in fj5l 

For this paper, we assume a concordance ACDM cos- 
mology, with Hf) = 70 km/s, JIa = 0.7, Vim = 0.3 and 
w = —1. For the given cluster rcdshift, the angular scale is 
8.46 kpc/". Magnitudes are given in the Vega system. 

2. Observations, data analysis, and results 

The cluster XMMU J0338.8-h0021 (hereafter X0338) was 
discovered within the XMM-A''ewio7T, Distant Cluster 
Project (XDCP), a serendipitous archival X-ray survey 
focussed on the identification of X-ray luminous systems 
at z > 0.8 (e.g. iMuUis et al.' 2005; B ohringer erall l2005t 



lFassbendeill2007t iFassbender et al.„20lTl ) 



2.1. X-ray selection with XMM-Newton 

The X-ray source associated with the cluster X0338 was ini- 
tially observed in 2002, in the XMM-Newton field with ob- 



servation identification number (OBSID) 0036540101 tar- 
geting the quasar SDSS033829.31-h002156 (z = 5.02) with 
a nominal exposure time of 22.9 ksec. It was detected as 
weak source with SAS v6.5 at an off-axis angle of 5.6', with 
a source significance of about 5 a and an extent significance 
of about 2.5 a. 

We have reprocessed the field with SAS vlO.O for a more 
accurate source characterization. After applying a strict 
two-step fiarc cleaning process for the removal of high back- 
ground periods, a clean net exposure time of 16.4/15.9 ksec 
remained for the two EMOS cameras and 7.8 ksec for the 
PN instrument. Figure[T] shows log-spaced X-ray contours 
(blue/cyan) of the cluster environment with significance 
levels of 2-16(7 derived from the adaptively smoothed com- 
bined images. For the flux measurement we applied the 
Towth curve analysis (GCA) method of iBohringer et al.l 
20001 ) in the soft 0.5-2 keV band and measured /x,500 — 
(3.0±1.8) X 10-15 ergs-icm-2 in a 40" (~i?5oo) aperture. 
At the cluster redshift, this translates into a 0.5-2 keV X- 
ray luminosity of ix,500 — (4.1 ± 2.4) x lO'*'^ erg/s. Owing 
to the faintness of the source at the limit of detectability, 
with ^33 net counts within an aperture of 24", the deter- 
mination of additional structural or spectral parameters is 
currently not feasible and the extended nature of the source 
is tent ative. However, ba sed on the M—Lx scaling rela- 
tion inlPratt et al.l (l2009f). w ith a modified redshift evolu- 
tion (jFassbender et alT ( 201lf ): Reichert et al., submitted), 
a flrst rough estimate of the expected temperature (r^^*) 
and mass (M500) can be obtained, with T_|?*^2.5keV and 
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-^'^500 ^ 0.8 X 10^'* Mq, corresponding to a total mass es- 
timate of M200 ~ 1-2 X IO^^Mq under the assumption of 
a NFW profile. According to the above Tl** value, a bolo- 
metric luminosity ix°500 — (1-1 iO. 6) x 10'*'* erg/s is finally 
estimated. We cannot rule out that a non-negligible part of 
the X-ray flux is due to contaminating point sources and 
thus the reported luminosities have to be taken as upper 
limits. All the estimated parameters based on the X-ray 
fiux are reported in Table[T] 



Table 1. Main properties of XMMU J0338.8+0021 



RA 




03:38:49.5 




DEC 




+00:21:08.1 




z 




1.490±0.009 




/x(0.5- 


- 2keV) 


(3.0 ±1.8) x 


10~'*ergs-'cm-^ 


Lx{0.5 


- 2keV) 


(4.1 ± 2.4) X 


10''=' erg s"^ 


r hoi 




(1.1 ±0.6) X 


10"*'' erg s"^ 


Mx,200 




~ 1.2 X 10^* 


Mq 


R200 




~ 590 kpc 


rpest 

J^x 




~ 2.5 keV 





2.2. Near-infrared follow-up imaging 

X0338 was photometrically confirmed as a significant over- 
density of very red galaxies coincident with the X-ray posi- 
tion using the wide-field (15.4' x 15.4') near-infrared (NIR) 
camera OMEGA2000 at the Calar Alto 3.5m telescope 
([Bailer- Jones et aO l2000l) . We obtained medium-deep H- 
band (50min) and z-band (53min) images of the X0338 
field on January 5, 2006 under good, but non-photometric, 
conditions with a seeing of 1.2". Re-observations with the 
z-filter were performed (5min) on October 30, 2006 in pho- 
tometric condi tions for the targ et and designated SDSS 
standard stars (jSmith et al.ll2002 ) in order to allow a pho- 
tometric calibration. The data wer e reduced with a desig- 
nated OME GA2000 NIR pipeline (lFassbendeij |2007l). The 
Sextractor (jBertin fc Arnoutjll996[ ) photometrv was cal- 
ibrated to the Ve ga system using 2MASS point sources 
(jCutri et al.ll2003f ) in H, and the SDSS standard star ob- 
servations in z, cross-calibrated with SDSS photometry in 
the science field. The limiting (Vega) magnitudes (50% 
completeness) were determined to be i?iim ~ 21.2 and 
ziini^23.1, corresponding to the expected apparent magni- 
tudes of passively evolving galaxies at Z'-^1.5 with absolute 



magnitude of M*+1.3 in H and M* in z. Figure[T] (left) 
shows an H-band image of the cluster field (4.5' x 4.5'), with 
overlaid X-ray contours (blue), density contours of color 
selected galaxies with 2.6 <z— H< 3.5 (red), and marked 
spectroscopic cluster members (small black circles) . A zoom 
onto the core region is shown in Fig.[Tl right. The z— H ver- 
sus H color magnitude diagram (CMD) of the field is shown 
in Fig.[2](top). 



3. Spectroscopic analysis 

3.1. Data reduction 

X0338 was observed with the F0RS2-VLT spectrograph 
on December 2009 (program ID:84.A-0844) with seeing 
conditions varying between 0.7"and 1.1". A single MXU 
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Fig. 2. Top: Color-magnitude diagram of the cluster field. 
All galaxies with Dproj < 30" are displayed in red. The 
seven secure spectroscopic cluster members are marked by 
black squares (with blue dots for Dproj > 30") and are all 
located within a color range 2.6<z— H<3.5 (green dot- 
ted lines), with an average color of z— H^3. The inter- 
lopers found in the above color cut are marked by black 
crosses. Two simple stellar population (SSP) models for 
the cluster redshift are overplotted at z— H~3, with solar 
metallicity and formation redshift of Zf = 5 (red dashed 
line) or Zf =3 (blue dashed). The vertical dotted blue 
line indicates the apparent magnitude of an L* galaxy at 
z = 1.49, black dashed lines show the 50%-completeness 
limit. Bottom: Composite rest frame spectrum of the seven 
members galaxies. The main features are marked by dot- 
ted lines and the total noise spectrum is shown in blue. 
An LRG template based o n z ~ 0.5 luminous red galax- 
ies (jEisenstein et al.l l2003l ) is overplotted in red to show 
the good match with the low-z template. The spectrum is 
smoothed with a 7 pixel boxcar filter. 
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Fig. 3. 1-D spectra of the seven confirmed members of XMMUJ0338. 8+0021, smoothed with a 7 pixel boxcar filter. For 
each spectrum, an identification label (ID) and the measured redsliift is reported. For the brightest component of the 
BCG (ID:1-A) the same LRG spectrum template used in Fig 12] (bottom) is overlaid in blue for comparison. In red, the 
sky features in absorption (top) and emission (bottom) arc shown. The position of the main spectral features, shifted 
according to the corresponding redshift values, are marked by the vertical dashed lines. 



Table 2. Properties of spectroscopic members of XMMU J0338. 8+0021 



ID 


RA 


DEC 


z 


0"z 


Projected Distance 


H 


z-H 


Main spectral 












(arcsec) 


(kpc) 


(Vega) 


(Vega) 


features 


1-A 


03:38:48.048 


+00:21:08.06 


1.5031 


0.0006 


20.8 


175.8 


18.25 


3.30 


Mgll, Fell, CaH/K 


5 


03:38:49.159 


+00:21:32.77 


1.5024 


0.0005 


24.9 


210.5 


20.19 


2.97 


Mgll, Fell, CaK 


7 


03:38:50.483 


+00:21:55.65 


1.4976 


0.0007 


49.8 


421.4 


19.17 


2.81 


Mgll, Fell (weak), [Oil] (weak) 


6 


03:38:48.449 


+00:21:45.78 


1.4899 


0.0006 


40.4 


341.9 


19.49 


3.54 


Fell, CaH/K 


4 


03:38:49.286 


+00:21:20.43 


1.4788 


0.0013 


12.6 


106.7 


19.53 


2.84 


CaH/K 


11 


03:38:52.357 


+00:22:43.76 


1.4736 


0.0005 


104.9 


887.7 


19.04 


2.94 


Mgll, [Oil] (weak), [Nelll], CaH/K 


14 


03:38:46.998 


+00:23:27.98 


1.4730 


0.0007 


144.5 


1222.2 


19.66 


3.33 


Mgll, [Oil] (weak), CaH/K 



slit-mask (field of view 6.8' x 6.8') with a total of 42 
1" width slits was used. The chosen instrument setup with 
the 3001 grism without blocking filter provides a wave- 
length coverage of 5800 - 10500A, with a resolution of 
R ~ 660. The total net integration time was 8.4ks. Raw 
data were reduced with a new F0RS2 adaption of the 
VIMOS Interac t ive P ipeline Graphical Interface (VIPGI, 
IScodeggio et al.l (|2005[ )^. which includes all the standard re- 
duction steps: bias subtraction, flat-field corrections, image 
stacking and wavelength calibration by means of a Helium- 
Argon reference spectrum (Nastasi et al., in prep.). 



A total of 45 1-dimensional spectra were obtained, with 
a final calibration uncertainty of ~lA. 

3.2. Redshift measurements 

In order to assess the redshift for each extracted spectrum, 
the graphical VIPGI tools were initially used for a quick vi- 
sual inspection. Then, the EZ software ([Garilli et al.ll20l"o[) 
and the RVSAO/IRAF package were also used with a set of 
different galaxy spectra templates. The analysis confirmed 
the existence of seven galaxies with concordant redshifts in 
the range 1.473<z< 1.503, six of which with a projected 
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Fig. 4. Left: Close-up z-band image of the BCG of XMMUJ0338.8-h0021, smoothed with a 0.7 arcsec gaussian filter. The 
contours indicate linearly spaced isophotes of the image cutout and the blue circle refers to an area with 3"of radius. The 
presence of a double-core system is clearly visible, revealing the brightest galaxy is experiencing an ongoing assembling 
process. Right: 2-D reduced spectra {top) and the 1-D extracted ones {bottom) for the two BCG components (ID: 1-A, 
1-B). In the bottom panel, the ID number and the measured redshift for each spectrum is reported in the top left corner, 
and the corresponding best-fit LRG template is overlaid in blue. Because of the low signal-to-noisc ratio obtained for 
ID: 1-B (whose spectrum is barely distinguishable in the top panel), its redshift value is only tentative. Both spectra are 
smoothed with a 7 pixel boxcar filter. 



distance from the X-ray emission peak Dproj < 2' and rest 
frame velocity offsets of |Az;| < 2000 km/s with respect to 
the median redshift. For the BCG, we extracted the spectra 
of both components ID:1-A and ID:1-B (Fig. El right), dis- 
tinguishable in the 0.86" seeing z-band preimage of Fig. 2] 
(left). For the second component (ID:1-B) only a tentative 
redshift measurement oi z = 1.5001 can be obtained ow- 
ing to its low signal-to-noise ratio. The extracted spectra 
and the main properties of the seven confirmed members 
are shown in Fig. [3] and reported in Table[2l respectively. 
Several spectral features are identifiable, albeit strong emis- 
sion and absorption sky lines affect the quality of the in- 
dividual spectra. This allows us to confirm the presence of 
a bound system of, at least, seven galaxies with a median 
redshift of z = 1.4 90 ib 0.009. Using the biweight estimator 
(|Beers et al.lll990[) we obtain a crude velocity dispersion 
estimate of cr^ « 670 km/s in the cluster rest frame, al- 
beit with an associated large uncertainty owing to the poor 
statistics. 

In order to enhance the faintest features, an average, 
noise-weighted composite spectrum of the seven members 
was created. The result is shown in Fig. [3] (bottom) and 
compared with an LRG template spectrum (red), derived 
for 726 luminous red galaxies at 0.47 < z < 0.55 by 
lEisenstein et all (|2003[ ). 



4. A young cluster in formation? 

From the CMD in Fig. [2] (top) we note that the red sequence 
of X0338 appears well populated. The 11 red galaxies, 
with colors 2.6 < z— H < 3.5 and brighter than the complete- 



ness limit, have a median (average) color of z— 11=2. 9 mag 
(3.0 mag), which is consistent with the model prediction of 
z— H~3.0 mag for a stellar formation redshift of Zf =3, 
which implies an evolved red sequence p opulation, as 
seen i n other distant c lusters (XMMUJ1229. ISantos et al.l 
(|2009D : XMMUJ2235, IStrazzullo et"aI1 (|2010D ). From the 
observed color spread of 5°J'_ffj « 0.25 mag, with an assumed 
flat slope and a Monte-Carlo estimated contribution of 
the photometric errors of (0.16±0.04) mag, we obtain a 
first order estimate of the intrinsic red sequence scatter of 
^T-H ^ (0.19±0.04) mag. This is significantly larger than in 
the above cases, e.g. 0.08 for XMMUJ22 35, but consistent 
with the findings of iHilton et"all (|2009D for their z=1.46 
cluster. The estimate for the above intrinsic scatter does 
not chan ge if a non-flat slope (e. g. the one of XMMUJ2235 
given bv IStrazzullo et ahl (|2010[ )) for the fit is used. This 
larger scatter can be explained with stellar formation age 
differences of the bright end of the red galaxy popul a tion o f 
the order of 0.5- 1 Gyr, as discussed in lHilton et al.l (|2009f ). 

The seven spectroscopic members show weak or no 
signs of ongoing star formation activity, with a mea- 
sured equivalent width of the [Oil] line in their spectra 
of |i^W^[07/]l < 5^- Assuming the above spectral feature as 
a proxy of their star formation ratcQ, the measured values 
for EW[oii] indicate a very moderate intensity for such ac- 
tivity. The existence of a real, albeit weak, [Oil] emission 
activity in the X0338 members is confirmed by the pres- 
ence of such a spectral feature in the composite spectrum 
of FiglH bottom. 

^ The presence of the [NeIII]A3870 line for the ID:11 galaxy, 
however, suggests the presence of an obscured AGN for it. 
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On the red sequence we find a magnitude gap be- 
tween the BCG and the second ran ked member o f 
Ami2 w 0.8 mag. From the study of .Smith et al.l (|2010[ ) 
in the local universe, this points towards a dynami- 
cally evolved cluster, while the apparent offset of the 
BCG with respect to the center of X-ray emission 
i^prof ~ 175 kpc) is a typical sign of a dynamically young 
clusters (jHaarsma et al.ll2010[ ). This is in line with the fact 
that the cluster shows indications of asymmetric X-ray 
emission with a low central surface brightness. A close-up 
look at the BCG, shown in Fig.|4] (left), shows a double- 
object morphology, indicating an active mass assembly pro- 
cess via a merger event (implying that the relatively large 
magnitude gap is a very recent feature) . Because of the lack 
of detectable emission lines in ID:1-A and 1-B spectra, we 
can state that this coalescence is likely a dry merger. The 
overall color and magnitude of the BCG are consistent with 
the ones found in red, passive, brightest galaxies residing in 
low-redshift clusters. The case of X0338 seems to suggest 
that we may be finally witnessing the active mass assembly 
epoch of BCGs, in contr ast to the evolved counterparts ob- 
serve d at z < 1.4 (e.g. IWhilev et all 120081: ICollins et all 
|2009() . Our observations support the hierarchical assem- 
bly scenario of BCGs, albeit at significantly earlier epoch 
than predicted from the c osmological simulation results of 
iDe Lucia fc BlaizotI (|2007D . 

The large scatter of the red sequence, the X-ray mor- 
phology and the spatial distribution of color-selected red 
galaxies point towards a younger cluster compared to e.g. 
XMMUJ2235, a scenario that will be critically tested by 
planned deeper observations. 

5. Summary and conclusion 

Wc summarize our results as follows: 

• We reported on the study of a newly discovered galaxy 
cluster XMMU JG338.8-I-0G21 at z = 1.490. This sys- 
tem was selected as an extended X-ray source serendip- 
itously detected within the XDCP. 

• By means of spectroscopic follow-up we can confirm 
seven cluster members, six of which with Dproj < 1 Mpc 
and a rest frame velocity offset of | Aw| < 2000 km/s from 
the median redshift. 

• The H-l-z band data reveal the presence of a well pop- 
ulated, albeit quite spread, red sequence in the CMD, 
with a median color of z— H = 2.9 mag and an intrin- 
sic color scatter of S'^J^'^ ^ 0.19. Comparisons with SSP 
models suggest that the stellar populations of red se- 
quence galaxies span the formation epoches in the red- 
shift range 3^2/^5. 

• A BCG is clearly identifiable on the red sequence and 
found in an active mass assembly phase, likely via a dry 
merging process. 

• Our data overall suggest that X0338 is a young galaxy 
cluster not yet fully evolved, as the faint X-ray emission 
appears to be morphologically disturbed and the large 
area over which its galaxies are found possibly indicates 
the presence of infalling structures. 

Currently, the present evidence for X0338 suggests 
partly an evolved cluster scenario and at the same time 
signs of a dynamically young system. More detailed stud- 
ies are needed for a more complete characterization of this 
intriguing system and its components. 
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